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Starting from a solution of the problem of a mechanical oscillator coupled to a scalar field inside 
Ph> a reflecting sphere of radius i?, we study the behaviour of the system in free space as the limit of 



Abstract 



an arbitrarily large radius in the confined solution. From a mathematical point of view we show 



5^ 

that this way of facing the problem is not equivalent to consider the system a priori embedded 
in infinite space. In particular, the matrix elements of the transformation turning the system to 
principal axis, do not tend to distributions in the limit of an arbitrarily large sphere as it should 



be the case if the two procedures were mathematically equivalent. Also, we introduce "dressed" 
coordinates which allow an exact description of the oscillator radiation process for any value of 
the coupling, strong or weak. In the case of weak coupling, we recover from our exact expressions 
the well known decay formulas from perturbation theory. 
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1 Introduction 



Since a long time ago the experimental and theoretical investigations on the polarization of atoms 
by optical pumping and the possibility of detecting changes in their polarization states has al- 
lowed the observation of resonant effects associated to the coupling of these atoms with strong 
radiofrequency fields ||]. As remarked in [Q], the theoretical understanding of these effects using 
perturbative methods requires the calculation of very high-order terms in perturbation theory, 
what makes the standard Feynman diagrams technique practically unreliable in those cases. The 
trials of treating non-perturbativelly such kind of systems consisting of an atom coupled to the 
electromagnetic field, have lead to the idea of " dressed atom" , introduced in refs M and || . This 
approach consists in quantizing the electromagnetic field and analyzing the whole system consist- 
ing of the atom coupled to the electromagnetic field. Along the years since then, this concept 
has been extensively used to investigate several situations involving the interaction of atoms and 
electromagnetic fields. For instance, atoms embedded in a strong radiofrequency field background 
in refs. || and M, atoms in intense resonant laser beans in ref. [0] or the study of photon correla- 
tions and quantum jumps. In this last situation, as showed in refs. H, H] and [TD|, the statistical 



properties of the random sequence of outcoming pulses can be analyzed by a broadband photode- 
tector and the dressed atom approach provides a convenient theoretical framework to perform this 
analysis. 

Besides the idea of dressed atom in itself, another aspect that desserves attention is the non- 
linear character of the problem involved in realistic situations, which implies, as noted above, in 



very hard mathematical problems to be dealt with. An way to circunvect these mathematical 
difficulties, is to assume that under certain conditions the coupled atom-electromagnetic field 
system may be approximated by the system composed of an harmonic oscillator coupled linearly 
to the field trough some effective coupling constant g. 

In this sense, in a slightly different context, recently a significative number of works has 
been spared to the study of cavity QED, in particular to the theoretical investigation of higher- 



generation Schrodinger cat-states in high-Q cavities, as has been done for instance in [II]. Linear 
approximations of this type have been applied along the last years in quantum optics to study 
decoherence, by assuming a linear coupling between a cavity harmonic mode and a thermal bath 
of oscillators at zero temperature, as it has been done in |L2] and ||13|| . To investigate decoherence 
of higher generation Schrodinger cat-states the cavity field reduced matrix for these states could 
be calculated either by evaluating the normal-ordering characteristic function, or by solving the 
evolution equation for the field-resevoir state using the normal mode expansion, generalizing the 



analysis of |12[ and [13 |. 

In this paper we adopt a general physicist's point of view, we do not intend to describe the 
specific features of a particular physical situation, instead we analyse a simplified linear version 
of the atom-field system and we try to extract the more detailed information we can from this 
model. We take a linear simplified model in order to try to have a clearer understanding of what 
we believe is one of the essential points, namely, the need of non-perturbative analytical treatments 
to coupled systems, which is the basic problem underlying the idea of dressed atom. Of course, 
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such an approach to a realistic non-linear system is an extremelly hard task and here we make 
what we think is a good agreement between physical reality and mathematical reliability, with the 
hope that in future work our approach could be transposed to more realistic situations. 

We consider a non relativistic system composed of a harmonic oscillator coupled linearly to 
a scalar field in ordinary Euclidean 3-dimensional space. We start from an analysis of the same 
system confined in a reflecting sphere of radius R, and we assume that the free space solution to 
the radiating oscillator should be obtained taking a radius arbitrarily large in the independent 
quantities. The limit of an arbitrarily large radius in the mathematics of the confined system 
is taken as a good description of the ordinary situation of the radiating oscillator in free space. 
We will see that this is not equivalent to the alternative continuous formulation in terms of 
distributions, which is the case when we consider a priori the system in unlimited space. The 
limiting procedure adopted here allows to avoid the inherent ambiguities present in the continuous 
formulation. From a physical point of view we give a non-perturbative treatment to the oscillator 
radiation introducing some coordinates that allow to divide the coupled system into two parts, the 
"dressed" oscillator and the field, what makes unecessary to work directly with the concepts of 
"bare" oscillator, field and interaction to study the radiation process. These are the main reasons 
why we study a simplified linear system instead of a more realistic model, to make evident some 
subtleties of the mathematics involved in the limiting process of taking a cavity arbitrarily large, 
and also to exhibit an exact solution valid for weak as well as for strong coupling. These aspects 
would be masked in the perturbative approach used to study non-linear couplings. 
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We start considering a harmonic oscillator q (t) of frequency u coupled linearly to a scalar 
field 0(r, t), the whole system being confined in a sphere of radius R centered at the oscillator 
position. The equations of motion are, 



q (0 + u 2 q (t) = 2n^gc I* dh^r, t)6(r) (1) 

Jo 

^" V20(r ' t) = 27r ^ go(t)5(r) (2) 
which, using a basis of spherical Bessel functions defined in the domain < |r| < R, may be written 
as a set of equations coupling the oscillator to the harmonic field modes, 

oo 

q (t) + u%q Q (t) =r]J2^iQi(t) (3) 

<ii(t) + u*qi(t) =r]Uiqo(t). (4) 

In the above equations, g is a coupling constant, r] = ^2gAu and Au = nc/R is the interval 
between two neighbouring field frequencies, u>i+i — Ui = Au = nc/R. 

2 The transformation to principal axis and the eigenfre- 
quencies spectrum 

2.1 - Coupled harmonic Oscillators 

Let us consider for a moment the problem of a harmonic oscillator go coupled to iV other 
oscillators. In the limit iV — > oo we recover our original situation of the coupling oscillator-field 



after redefinition of divergent quantities, in a manner analogous as renormalization is done in field 
theories. In terms of the cutoff iV the coupled equations (§) and @ are simply rewritten taking 
the upper limit N instead of oo for the summation in the right hand side of Eq. ([3D and the system 
of N + 1 coupled oscillators go {Qi} corresponds to the Hamiltonian, 
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N 



Po + ^oVo + Yl Pk + ufal - ZvUkqoQk 



k=l 



(5) 



The Hamiltonian (0) can be turned to principal axis by means of a point tranformation, 



(6) 



performed by an orthonormal matrix T = (tV), fi = (0, k), k = 1,2, ...N, r = 0,...N. The 
subscript and k refer respectively to the oscillator and the harmonic modes of the field and r 
refers to the normal modes. The transformed Hamiltonian in principal axis is 



N 



z r=0 



(7) 



where the fi r 's are the normal frequencies corresponding to the possible collective oscillation modes 
of the coupled system. 

Using the coordinate transformation = t r ^Q r in the equations of motion and explicitly 



making use of the normalization condition J2n=o(t r n) 2 = 1, we get 



V^k r 
'0- 



(8) 



'o 



1 + 



N 



(9) 



and 

N 2 

-o 2 -^ 2 £-^V (io) 

k=l U k-^ L v 

There are N + 1 solutions Q r to Eg . (|T0|) , corresponding to the N + 1 normal collective oscillation 
modes. To have some insight into these solutions, we take fl r = Q in Eq.([H]) and transform the 
right hand term. After some manipulations we obtain 

k=l u k - " 

It is easily seen that if uo 2 > Nrf Eq.([TT|) yelds only positive solutions for Q 2 , what means that the 
system oscillates harmonically in all its modes. Indeed, in this case the left hand term of Eq. QTTj ) 
is positive for negative values of Q 2 . Conversely the right hand term is negative for those values 
of Vt 2 . Thus there is no negative solution of that equation when > Nrj 2 . On the other hand it 
can be shown that if uo 2 < Nrj 2 , Eq.(|ll]) has a single negative solution Vt 2 _. In order to prove it let 
us define the function 

m 2 ) = k 2 - n v 2 ) -n 2 - v 2 J2 (12) 

k=i u k - * l 

Accordingly Eq . (|Tl|) can be rewritten as I(fi 2 ) = 0. It can be noticed that I(Q 2 ) — >• oo as 
Q 2 — > — oo and 

I(Q 2 = 0) = u 2 - Nrf < (13) 

Furthermore /(^ 2 ) is a monotonically decreasing function in that interval. Consequently I(Q 2 ) = 
has a single negative solution when uo 2 < Nrf as we have pointed out. This means that there 
is an oscillation mode whose amplitude varies exponentially and that does not allows stationary 
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configurations. We will not care about this last situation. Thus we assume > Nrj 2 and define 



the renormalized oscillator frequency uj [14 



oo = sfu% - N V 2 . (14) 



In terms of the renormalized frequency Eq.(|Tl]) becomes, 

A n 2 



-2 



"--V = v 2 Y.-r^- (15) 



/, =1 U k - --r 



From Eqs. (|8]), (|9]) and (|15D, a straightforward calculation shows the orthonormality relations for 
the transformation matrix (t r ). 

We get the transformation matrix elements for the oscillator-field system by taking the limit 
N — > oo in the above equations. Recalling the definition of r\ from Eqs. (|3]) and ([|), we obtain 
after some algebraic manipulations, from Eqs. (|i5|) , (|8]) and (^|), the matrix elements in the limit 
N -> oo, 

tl = ^ (16) 



and 



U! 



2.2 - The eigenfrequencies spectrum 

Let us now return to the coupling oscillator-field by taking the limit iV — >• oo in the relations of 
the preceeding subsection. In this limit it becomes clear the need for the frequency renormalization 
in Eq . fllil) . It is exactly the analogous of a mass renormalization in field theory, the infinite ujq is 



chosen in such a way as to make the renormalized frequency uo finite. Remembering Eq.(^) the 
solutions with respect to the variable Q of the equation 

w - fi =^r£^' (18 > 

give the collective modes frequencies. We remember Uk = k j^, k — 1,2, and take a positive x 
such that fl = Xj^. Then using the identity, 

00 x 2 1 



Eq . (|P8|) may be rewritten in the form, 

c 1 .„ ito 2 . 

cotgirx = —x-\ (1 ). (20) 

Kg tix ngc 

The secant curve corresponding to the right hand side of the above equation cuts only once each 
branch of the cotangent in the left hand side. Thus we may label the solutions x r as x r — r + e r , 
< e r < 1, r = 0,1, 2..., and the collective eigenfrequencies are, 

n r = (r + e r )— , (21) 

the e's satisfying the equation, 

Vtl — uo 2 c . . 

co «^> = -fer + tw- (22) 

The field </>(r, t) can be expressed in terms of the normal modes. We start from its expansion 
in terms of spherical Bessel functions, 

00 

0(r,t)=c5> fc (t)0 fc (r), (23) 
fe=i 
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where 



sin— r| 



Mr) = —ff=- (24) 
Using the principal axis transformation matrix together with the equations of motion we obtain 
an expansion for the field in terms of an orthonormal basis associated to the collective normal 
modes, 

oo 

0(r,t) = c^Q s (t)$ s (r), (25) 

s=0 

where the normal collective Fourier modes 

*.w=£n f ^ (26) 

k rv27rit 

satisfy the equation 

(-5* - A)0 s (r) = 27r^(rK (27) 



which has a solution of the form 



^Vfdk^i-^ (28) 



To determine the phase S s we expand the right hand term of Eq.(^) and compare with the formal 
expansion (pf). This imply the condition 

sin(—R - 5 S ) = 0. (29) 

c 

Remembering from Eq.(pT|) that there is < e s < 1 such that fl s = (s + e s )j^, it is easy to show 
from the condition in Eq. (|27|) that the phase < 5 S < tt has the form 

S s = e s 7r- (30) 
10 



Comparing Eqs.(|24|) and (^) and using the explicit form ([16]) of the matrix element we obtain 
the expansion for the field in terms of the normal collective modes, 

A ) ( Yt ) = -^y Q s sin{^\v\-5 S ) 

3 The limit R —> oo - mathematical aspects 

3.1 - Discussion of the mathematical problem 

Unless explicitly stated, in the remaining of this paper the symbol R — > oo is to be understood 
as the situation of a cavity of fixed, arbitrarily large radius. In order to compare the behaviour 
of the system in a very large cavity to that it would be in free space, let us firstly consider the 
system embedded in an a priori infinite Euclidean space; in this case to compute the quantities 
describing the system means essentially to replace by integrals the discrete sums appearing in the 
confined problem, taking direcltly R = oo. An alternative procedure is to compute the quantities 
describing the system confined in a sphere of radius R and take the limit R — > oo afterwards. This 
last approach to describe the system in free space should keep in some way the "memory" of the 
confined system. To be physically equivalent one should expect that the two approachs give the 
same results. We will see that at least from a mathematical point of view this is not exactly the 
case. We remark that a solution to the problem of a system composed of an oscillator coupled to 



a field in free space, is already known since a long time ago |15] in the context of Bownian motion. 



This solution is quite different from ours, in the sense that it not concerns the system confined to 
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a box and also that it is limited to the dipole term from the multipolar expansion to the field. 

In the continuous formalism of free space the field normal modes Fourier components (analogous 
to the components S in Eq.fl2"6|)) are, 

, , r°° oo sin-\r\ , 

(Pn = h(Q) / duo . 09 , c ,' (32) 



o 



uo 2 -il 2 Irl 



h(Q) = 2 £ = (33) 

'(Q 2 - uo 2 ) 2 + 7cg 2 fl 2 



where 



and where the we have taken the appropriate continuous form of Eqs.fllBD and QI7|). Splitting 
ujjiuj 2 — Q 2 ) into partial fractions we get 

<f> n = h(Sl) / duo (34) 



CO 



oo — Q Irl 



The pole at lo = f2 prevents the existence of the integral in Eq. (|34]) . The usual way to 
circumvect this difficulty is to replace the integral by one of the quantities, 

f+°° 1 sm-|r| f+°° sin- Irl 

Lim e ^Q I duo — — - — — — n — = / duoo±{uo — \l) — j-^ — , (35) 



oc 



uo — (Q ± ie) | x* | 



where 



5 ± {uo -Q) = -P{—^—)±i5{uo -Q), (36) 

7T UO — \l 

with P standing for principal value. In our case this redefinition of the normal modes Fourier 
components may be justified by the fact that both integrals in Eq. ([35]) are solutions of the equations 
of motion ([[]) and (Q) for r / 0, and so the solution should be a linear combination of them. The 
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situation is different if we adopt the point of view of taking the limit R — > oo in the solution of the 
confined problem. In this case the Fourier component <ftn is obtained by taking the limit R —>■ oo 
in the expression for the field, Eq(p8|), what allows to obtain an uniquely defined expression to the 
normal modes Fourier components, to each <pn corresponding a phase 5n (the limit R — > oo of 5 S 
in Eq-dUD given by 

cotdn = . (37) 

ng si 

Also, comparing Eqs.(|35|), (|3~6|) and (^) we see that the adoption of the continuous formalism is 
equivalent to assume that in the limit R — ■> oo the elements if of the transformation matrix should 
be replaced by 5 + (uj — Q) or by 5-(uj — Q). This procedure is, from a mathematical point of view, 
perfectly justified but at the price of loosing uniqueness in the definition of the field components. 

If we take the solution of the confined problem and we compute the matrix elements if for R 
arbitrarily large, we will see in subsection 3.2 that these elements do not tend to distributions in 
this limit. As R becomes larger and larger the set of non- vanishing elements if concentrate for 
each i in a small neighbourhood of uji. In the limit R — > oo the whole set of the matrix elements 



tf contains an arbitrarily large number of elements quadratically summables [0 . For the matrix 
elements we obtain a quadratically integrable expression. 

In the continuous formulation the unit matrix, corresponding to the absence of coupling, has 
elements E„ = 5(u — Q), while if we start from the confined situation, it can be verified that in the 
limit g — > 0, R —>■ oo, the matrix T = (t s ) tends to the usual unit matrix of elements E u> q = 5 W) q. 

The basic quantity describing the system, the transformation matrix T = (t s ) has, as we 
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will see, different properties in free space, if we use the continuous formalism or if we adopt the 
procedure of taking the limit R — > oo from the matrix elements in the confined problem . In 
the first case we must define the matrix elements v£ linking free field modes to normal modes, 
as distributions. On the other side adopting the second procedure we will find that the limiting 
matrix elements Lim^^ t\ are not distributions, but well defined finite quantities. The two 
procedures are not equivalent, the limit R — > oo does not commute with other operations. In this 
note we take as physically meaningfull the second procedure, we solve first the problem in the 
confined case (finite R) and take afterwards the limit of infinite (in the sense of arbitrarily large) 
radius of the cavity. In the next subsection we perform a detailed analysis of the limit R — > oo of 
the transformation matrix (t r ). 

3.2 - The transformation matrix in the limit R — > oo 

From Eqs. ( pp and ( |TTD we obtain for R arbitrarily large, 



where we have used the fact that in this limit Aoo = Afi = The matrix elements tn are 
quadratically integrable to one, J(t2) 2 dfl = 1, as may be seen using Cauchy theorem. 

For R arbitrarily large (Aa> = ^ —>■ 0), the only nonvanishing matrix elements t\ are those 
for which Ui — Q r ~ Au. To get explicit formulas for these matrix elements in the limit R — > oo 




(38) 



and 



2guJkAuj 



(39) 
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let us consider R large enough such that we may take Auj ~ AQ and consider the points of the 
spectrum of eigenfrequencies Q inside and outside a neighbourhood rj (defined in Eqs.@ and 
of uJi. We note that R > — implies | > Au, then we may consider R such that the right (left) 
neighbourhood ^ of Ui contains an integer number, k, of frequencies Q r , 



^ = 2=^. (40) 

If R is arbitrarily large we see from ( fTO| ) that | is arbitrarily small, but n grows at the same rate, 
what means firstly that the difference uji — Q r for the Q r J s outside the neighbourhood rj of uoi is 
abitrarily larger than Au, implying that the corresponding matrix elements t\ tend to zero (see 
Eq.(PD|)). Secondly all frequencies fl r inside the neighbourhood 77 of uJi are arbitrarily close to u>i, 
being in arbitrarily large number. Only the matrix elements t\ corresponding to these frequencies 
Q r inside the neighbourhood rj of uo,i are different from zero. For these we make the change of 
labels, 

r = i — n [uJi < f2 r < cjj) ; r = i + n (&i > £l r > u){ H — ), (41) 



a = — 9 r: - i (42) 



1,2,.... We get, from Eq.(g| 

t- = - 

'{p 2 - u 2 ) 2 + ir 2 g 2 ut 



and 



,i±n Tg^i 1 



±n = (43) 
(tt 2 -u 2 ) 2 + n 2 g 2 LJ 2n ±^ 



where e« satisfies Eq.(|22|) in this case, 



u 2 -u 2 



coi(7re .) = _J . (44) 
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Using the formula 



7r 2 cosec 2 (iT€j 



-* n=l 



[n + 6; 



+ 



'n - €if 



(45) 



it is easy to show the normalization condition for the matrix elements (f42|) and ( [43]) 



(*i) 2 + E(* 

and also the orthogonality relation, 



i— n\2 I /j.i+n\2 



(46) 



n=l 



Et;tj; = o(<^fc) 

r 

in the limit i? — > oo. 

3.3 - The transformation matrix in the limit g = 

From Eq. ([TBI) we get for arbitrary R, 

1, if fi r = cD; 
0, otherwise. 

From Eqs. (|42|) and ([43]) we see that the matrix elements t\ for % ^ r all vanish for g 
using Eqs.(^l|) we obtain for small g, 

2gQiUi 1 



(47) 



Lim g -> t r 



(48) 
0. Also, 



t 



or, expanding for small g from Eq. 



(n 2 - ^)(ui + tit) ej 



*K<? = 0) = 1 



(49) 



(50) 



We see from the above expressions that in the limit R — > oo the matrix (tp remains an 
orthonormal matrix in the usual sense as for finite i?. With the choice of the procedure of taking 
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the limit R — > oo from the confined solution, the matrix elements do not tend to distributions in 
the free space limit as it would be the case using the continuous formalism. All non- vanishing 
matrix elements t\ are concentrated inside a neighbourhood rj of u>i, their set is a quadratically 
summable enumerable set. The elements (tj) tend to a quadratically integrable expression. 

4 The Radiation Process 

We start this section defining some coordinates q' , q\ associated to the "dressed" mechanical 
oscillator and to the field. These coordinates will reveal themselves to be suitable to give an 
appealling non-perturbative description of the oscillator-field system. The general conditions that 
such coordinates must satisfy, taking into account that the system is rigorously described by the 
collective normal coordinates modes Q r , are the following: 

- In reason of the linear character of our problem the coordinates q' , q\ should be linear functions 
of the collective coordinates Q r 

- They should allow to construct ortogonal configurations corresponding to the separation of 
the system into two parts, the dressed oscillator and the field. 

- The set of these configurations should contain the ground state, r . 

The last of the above conditions restricts the transformation between the coordinates q' , fi = 
0,i = 1,2, ... and the collective ones Q r to those leaving invariant the quadratic form, 




(51) 



r 
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Our configurations will behave in a first approximation as independent states, but they will evolve 
as the time goes on, as if transitions among them were being in progress, while the basic configu- 
ration T represents a rigorous eigenstate of the system and does not change with time. The new 
coordinates q' describe dressed configurations of the oscillator and field quanta. 
4.1 - The dressed coordinates q'^ 

The eigenstates of our system are represented by the normalized eigenfunctions, 



nonw2- 



XQ,t) = U 



N n H n J^Q s ) 



r e -i ^ nsnst , 



where H Us is the n s -th Hermite polynomial, N Hs is a normalization coefficient, 



(52) 



N„ 



(53) 



and T is a normalized representation of the ground state, 



T = exp 



~v 2h 4 n ^h 



(54) 



To describe the radiation process, having as initial condition that only the mechanical oscillator, 
q be excited, the usual procedure is to consider the interaction term in the Hamiltonian written 
in terms of go, Qi as a perturbation, which induces transitions among the eigenstates of the free 
Hamiltonian. In this way it is possible to treat approximately the problem having as initial 
condition that only the bare oscillator be excited. But as is well known this initial condition is 
physically not consistent due to the divergence of the bare oscillator frequency if there is interaction 
with the field. The traditional way to circumvect this difficulty is by the renormalization procedure, 
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introducing perturbativelly order by order corrections to the oscillator frequency. Here we adopt 
an alternative procedure, we do not make explicit use of the concepts of interacting bare oscillator 
and field, described by the coordinates go and {qi}, we introduce "dressed" coordinates q' Q and 
{q'j} for, respectively the "dressed" oscillator and the field, defined by, 

/fi-ix/fa, (55) 



valid for arbitrary R, which satisfy the condition to leave invariant the quadratic form ( pT|) and 
where = Q, {ui}. In terms of the bare coordinates the dressed coordinates are expressed as, 

i = (XiwQv, ( 56 ) 

V 

where 



1 



a 



r 



As R becomes larger and larger we get for the various coefficients a in Eq. (CT) : 
a) from Eq . (|38|) . 

1 f°° 2gQ 2 y/UdQ 
LirriR^oo a 00 = —t= / j-rr r^- — Yl^a = A oo(u,g). (58) 

y/U Jo {\l z — UJ Z ) Z + 7T z g z ll z 

b) To evaluate and aoi in the limit R —>■ oo, we remember from the discussion in subsection 
3.2 that in the the limit R —> oo, for each i the only no n- vanishing matrix elements t\ are those 
for which the corresponding eigenfrequencies Q r are arbitrarily near the field frequency u>i. We 
obtain from Eqs. (153), (PE3) and 



ful (wf - u z y + 7r z g z uf ^ n 2 - ef - 
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L»mB_»oo a i = Lim^^Q — 5-^-^(2^ -5 2 ) ( 60 ) 

v4J (utf - ^j 2 ) 2 + Tx 2 g 2 uf ^ n 2 - ef €i 



and 



c) Since in the limit R —>■ 00 the only non-zero matrix elements t\ corresponds to Q r = cjj, the 
product vanishes for uj,i 7^ cj^. Then we obtain from Eqs. (|57D and (^) 

Limji^oo a ik = 5 ik . (61) 

Thus, from Eqs. (|56|) , (|6T|) , (|59|) , (|60|) and (|58|) we can express the dressed coordinates g^ in terms 
of the bare ones, g M in the limit R — > 00, 

go = A 00 g , (62) 

Qi = H- (63) 

It is interesting to compare Eqs.(^6|) with Eqs. (|62|) , (|63"D. In the case of Eqs. (|56"D for finite 
R, the coordinates q' and {q^} are all dressed, in the sense that they are all collective, both the 
field modes and the mechanical oscillator can not be separeted in this language. In the limit 
R — ► 00, Eqs.(|62"D and (|63| ) tells us that the coordinate q' describes the mechanical oscillator 
modified by the presence of the field in a indissoluble way, the mechanical oscillator is always 
dressed by the field. On the other side, the dressed harmonic modes of the field, described by the 
coordinates q[ are identical to the bare field modes, in other words, the field keeps in the limit 
R — ► 00 its proper identity, while the mechanical oscillator is always accompanied by a cloud of 
field quanta. Therefore we identify the coordinate q' as the coordinate describing the mechanical 
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oscillator dressed by its proper field, being the whole system divided into dressed oscillator and 
field, without appeal to the concept of interaction between them, the interaction being absorbed 
in the dressing cloud of the oscillator. In the next subsections we use the dressed coordinates to 
describe the radiation process. 

4.2 - Dressed configurations and the radiation process 

Let us define for a fixed instant the complete orthonormal set of functions, 

= n 

where q'^ = q' Q , q[, uj^ = u>, Ui and N K/1 and T are as in Eq. (|52]) . Using Eq.flolf) the functions 
(|64|) can be expressed in terms of the normal coordinates Q r . But since (|52|) is a complete 
set of orthonormal functions, the functions (Q) may be written as linear combinations of the 
eigenfunctions of the coupled system (we take t = for the moment), 

Vw.(?0= E ^T 1 1 ::(°)Kn,n 2 ...(Q,o), (65) 

nom — 

where the coefficients are given by, 

t: o °Z:;(o) = J dQVw.</W 2 ..., (66) 

the integral extending over the whole Q-space. 

We consider the particular configuration ip in which only one dressed oscillator q'^ is in its iV-th 
excited state, 

^o...n M o...W) = N^i^q'^To. (67) 
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r , 



(64) 



The coefficients fl6"6|) can be calculated in this case using Eqs.(|66"]), ( |64D and (^) with the help of 
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ml 



E(^) 2 



^( /y; )= E — , - ^o(t yQo)g mi (^Qi)... 



(68) 



We get, 
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T o"°.;W.. = (^T-T-) 5 (^) n °(^) ni - ; (69) 

I0Q.I01 



where the subscripts // = 0, i refer respectively to the dressed mechanical oscillator and the 
harmonic modes of the field and the quantum numbers are submited to the constraint rio+Tii+... = 
N. 

In the following we study the behaviour of the system with the initial condition that only the 
dressed mechanical oscillator q' be in the N-th excited state. We will study in detail the particular 
cases N — 1 and N = 2, which will be enough to have a clear understanding of our approach. 

- N — 1: Let us call the configuration in which only the dressed oscillator q'^ is in the first 
excited level. The initial configuration in which the dressed mechanical oscillator is in the first 
excited level is r°. We have from Eq. (|67"D , (|65|) flBTf ) and (j55|) the following expression for the time 
evolution of the first-level excited dressed oscillator q'^, 

r^ = J2r(tn(o), (70) 

V 

where the coefficients f^ u {t) are given by 

r^) = E^ e ~^ ( 71 ) 

s 
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That is, the initially excited dressed oscillator naturally distributes its energy among itself and 
all others dressed oscillators, as time goes on. If the mechanical dressed oscillator is in its first 
excited state at t = 0, its decay rate may evaluated from its time evolution equation, 

rS = £/ 0,, (*)ry(o). (72) 

V 

In Eq.(|72"D the coefficients f 0v {t) have a simple interpretation: remembering Eqs.(^) and (|63D, 
f 00 (t) and f 0l (t) are respectively the probability amplitudes that at time t the dressed mechanical 
oscillator still be excited or have radiated a field quantum of frequency cjj. We see that this for- 
malism allows a quite natural description of the radiation process as a simple exact time evolution 
of the system. Let us for instance evaluate the oscillator decay probability in this language. From 
Eqs.(|38[) and ([fl]) we get 

00 _ too ggov^ <m 

J W-y (W-U*)*+TT*gW W 

The above integral can be evaluated by Cauchy theorem. For large t (t » 4), but arbitrary 
coupling g, we obtain for the oscillator decay probability, the result, 

, r (t)| . = + ^) + e -^(^ + J ^ t) + ig!, (74) 

where a) = J Co 2 — In the above expression the approximation t » -r plays a role only in 

the two last terms, due to the difficulties to evaluate exactly the integral in Eq. ([73|) along the 
imaginary axis. The first term comes from the residue at Q = u + i^f and would be the same if 
we have done an exact calculation. If we consider the case of weak coupling, g << u, we obtain 
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the well known perturbative exponential decay law for the harmonic oscillator ||18||, 

|/ 00 (t)| 2 ^e-^\ (75) 

but we emphasize that Eq.([74|) is valid for all values of the coupling constant g, even large, it is 
an expression valid for weak as well as strong couplings. 
- N = 2 

Let us call T±[ the configuration in which the dressed oscillators q'^ and q' v are at their first 
excited level and the configuration in which q'^ is at its second excited level. Taking as initial 
condition that the dressed mechanical oscillator be at the second excited level, the time evolution 
of the state may be obtained in an analogous way as in the preceeding case, 

r§(t) = E \rit)f + 4 E f°" (t)f* (76) 

where the coefficients j m and f 0fl are given by ([7ip. Then it easy to get the following probabilities: 
Probability that the dressed oscillator still be excited at time t: 

Po(t) = \f°(t)\\ (77) 

probability that the dressed oscillator have decayed at time t to the first level by emission of a 
field quantum: 

p 1 (t) = 2ir°(t)i 2 (i-ir°(t)i 2 ) (78) 

and probability that the dressed oscillator have decayed at time t to the ground state: 

P 2 (t) = l-2|/ 00 (t)| 2 + |/ 00 (t)| 4 . (79) 
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Replacing Eq.(|T3|) in the above expressions we get expressions for the probabilities decays valid 
for any value of the coupling constant. In the particular case of weak coupling we obtain the well 
known perturbative formulas for the oscillator decay . 



P (t) « e~ 2n9 \ (80) 

Pi (t) « 2e -ff * t (l - e -,rat ) (81) 

and 

P 2 (t) « 1 - 2e" Kgi + e- 27Tgt . (82) 



5 Concluding Remarks 

In this paper we have analysed a symplified version of an atom-electromagnetic field system and 
we have tried to give the more exact and rigorous treatment we could to the problem. We have 
adopted a general physicist' s point of view, in the sense that we have rennounced to approach 
very closely to the real behaviour of a complicated non-linear system, to study instead a simple 
linear model. As a counterpart, an exact solution has been possible. Our dressed coordinates give 
a description of the behaviour of the system that is exact and valid for weak as well as for strong 
coupling. If the coupling between the mechanical oscillator and the field is weak, we recover the 
well known behaviour from perturbation theory. 
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6 In Memoriam 



This paper evolved from umpublished work we have done and discussions we have had, with Prof. 
Guido Beck when two of us (A.P.C.M. and N.P.A.) were his students at Instituto de Fisica Balseiro 
in Bariloche (Argentina), in the late sixties and the early seventies. We dedicate this article to his 
memory. 
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